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The interaction of optical and mechanical modes in nanoscale optomechanical systems has been widely 
studied for applications ranging from sensing to quantum information science. Here, we develop a plat¬ 
form for cavity optomechanical circuits in which localized and interacting 1550 11 m photons and 2.4 GHz 
phonons are combined with photonic and phononic waveguides. Working in GaAs facilitates manipula¬ 
tion of the localized mechanical mode either with a radio frequency field through the piezo-electric effect, 
or optically through the strong photoelastic effect. We use this to demonstrate a novel acoustic wave 
interference effect, analogous to coherent population trapping in atomic systems, in which the coherent 
mechanical motion induced by the electrical drive can be completely cancelled out by the optically-driven 
motion. The ability to manipulate cavity optomechanical systems with equal facility through either pho¬ 
tonic or phononic channels enables new device and system architectures for signal transduction between 
the optical, electrical, and mechanical domains. 


The interaction of optical and mechanical degrees of free¬ 
dom in chip-based, nanoscale systems has been studied in 
many contexts, ranging from cavity optomechanical sys¬ 
tems- 3 in which localized optical and mechanical modes are 
coupled, to waveguiding geometries in which propagating 
photons and phonons interact, as in stimulated Brillouin scat¬ 
tering' 4 '® (the latter has also been studied in travelling-wave 
microresonators 7 9 ). Recently there has been interest in com¬ 
bining the exquisite motion sensitivity of cavity optomechani¬ 
cal systems with the radio frequency (RF) signal processing 
functionality of electromechanical systems to enable wave¬ 
length conversion between microwave and optical domain^- 1 ' 17 
electrostatically-actuated optomechanical cavities for sensing 
and feedback cooling 11 2T anc | piezo-optomechanical cavi¬ 
ties 13 14 in which mechanical motion is driven by RF fields, 
with sensitive readout and coherent interference effects ob¬ 
servable in the optical domain. 

In this work, we develop a platform that combines the 
strong interaction achievable in a nanoscale cavity optome¬ 
chanical system with coherent control of the localized me¬ 
chanical mode along both the optical channel 1 ? ®and through 
coupling to propagating phonons using acoustic waveguides. 
Working in the GaAs material system allows us to combine 
its photoelastic effect, which strongly couples localized opti¬ 
cal and mechanical modes 17 SI, with its piezoelectric behav¬ 
ior, which enables the generation of surface acoustic waves 
from an RF drive 19 . We develop phononic crystal waveguides 
that guide acoustic waves that couple to the 2.4 GHz breathing 
mode of nanobeam optomechanical crystal cavities 18 20 21 and 
are sourced by RF-driven interdigitated transducers (IDT). 
Phase-sensitive optical measurements confirm the ability to 
initialize the mechanical resonator with arbitrary amplitude 
and phase, and to detect an average coherent intracavity 
phonon population of less than one. Finally, we take advan¬ 
tage of an optomechanical coupling rate go/2n « 1.1 MHz 
that is one to two orders of magnitude larger than that achieved 
in other piezo-optomechanical systems l3 -lSlto demonstrate an 


acoustic interference effect analogous to atomic coherent pop¬ 
ulation trapping, in which the RF-driven coherent motion of 
the mechanical resonator is cancelled by the optically-induced 
motion, and vice versa. Coherent manipulation of the me¬ 
chanical resonator through both optical and electrical chan¬ 
nels with equal facility, combined with on-chip routing pro¬ 
vided by phononic crystal waveguides, enables new applica¬ 
tions of integrated electro-optomechanical systems in signal 
transduction and sensing. 

Coupling Propagating and Localized Phononic Modes 

Figure |TJa) shows an overview of our optomechanical cir¬ 
cuit platfomPS. The system uses an inter-digitated transducer 
(IDT) to convert an applied RF drive to a propagating surface 
acoustic wave via the piezoelectric effect. This wave is routed 
using a phononic crystal waveguide and is butt-coupled to a 
nanobeam optomechanical crystal cavity (Fig. [ljc)-(d)). The 
optical interface to the cavity is a fiber taper waveguide, which 
is used to inject and extract light from the device. Finally, a 
second IDT can electrically detect surface acoustic waves out- 
coupled from the nanobeam cavity, or act as a second source 
of surface acoustic waves that couple to the cavity via a sec¬ 
ond phononic crystal waveguide. 

Figure |T] b) shows the optomechanical crystal cavity de¬ 
sign, which consists of an array of elliptical air holes in a 
suspended nanobeanPD. The quadratic grade of the lattice 
constant and major and minor axis diameters of the elliptical 
holes co-localizes 1550 nm photons and 2.4 GHz phonons to 
a micrometer length scale. The optomechanical coupling rate 
go/2n, which quantifies optical cavity frequency shift due to 
zero-point motion, is «1.1 MHz and is predominantly due to 
the photoelastic effect. This is amongst the highest reported 
for a cavity optomechanical system, and enables strong dy¬ 
namic back-action so that the optical field can coherently ma¬ 
nipulate the localized mechanical mode. For coherent control 
of the localized mechanical breathing mode in the RF domain, 
we need mechanisms to convert the RF voltage into a propa¬ 
gating acoustic wave and a waveguide geometry to route these 
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FIG. 1: Piezo-optomechanical circuits, (a) Optical 
micrograph of the device layout. An array of op¬ 
tomechanical cavities, each coupled to input and out¬ 
put phononic crystal waveguides, is placed between 
inter-digitated transducers (IDTs) that excite and/or 
detect surface acoustic waves. Optical coupling is 
done using a fiber taper waveguide, shown schemat¬ 
ically in gray, (b) Optomechanical crystal cavity de¬ 
sign. (Top) Variation of the lattice constant, major 
axis diameter, and minor axis diameter of the elliptical 
holes along the nanobeam; (Middle) Simulated nor¬ 
malized electric field amplitude for the 1550 nm op¬ 
tical mode. (Bottom) Simulated normalized displace¬ 
ment amplitude for the 2.4 GHz mechanical breathing 
mode, (c) Scanning electron microscope (SEM) im¬ 
age of a nanobeam optomechanical crystal cavity, in¬ 
cluding the adjacent input and output phononic crystal 
waveguides, (d) SEM image showing the transition 
from the input IDT (shaded in blue) to the phononic 
crystal waveguide and nanobeam cavity. 


waves and couple to the nanobeam cavity. We begin by con¬ 
sidering the acoustic wave generation. 

Surface acoustic waves (SAWs) have long been used for 
applications in signal processing and communication 23 24 . Al¬ 
though they exist on the surface of every material, piezoelec¬ 
tric materials are used most extensively because of the ease 
of SAW excitation using IDTs, which are interleaved metal 
fingers biased with opposite polarity 22 . To take into account 
the multilayer metal electrode and underlying epitaxial layer 
structure (GaAs on Alo. 7 Gao. 3 As), the IDT resonance fre¬ 
quency is determined from finite element method (FEM) sim¬ 
ulations. Figure[2ja) shows the simulated SAW cross-section. 

To route acoustic energy to the nanobeam optomechanical 
cavity, we use phononic crystal waveguides. The phononic 
crystal design is adapted from the periodic cross-structure- 23 , 
which supports a complete bandgap for acoustic waves 22 . The 
waveguide is a line defect geometry created by removing one 
row of crosses, resulting in a structure that supports laterally- 
confined, propagating acoustic modes. Figure [2| a) shows an 
FEM simulation of an acoustic wave excited by an IDT and 
propagating through the waveguide. The phononic crystal 
confines the mechanical energy to the line defect region, with 
little displacement beyond a couple of periods transverse to 
the propagation direction. 

We use FEM simulations to understand the coupling be¬ 
tween the localized and propagating mechanical modes. Fig¬ 
ure [ 2 ]b) shows the results of a simulation in which the waveg¬ 
uide mode at 2.25 GHz is launched into the nanobeam cav¬ 
ity, with excitation of the resonant mechanical breathing mode 
clearly observed. This coupling can be understood quantita¬ 
tively from the cavity’s acoustic reflection spectrum 22 . The 
finite intrinsic Q m of the cavity in fabricated devices is mod¬ 
eled by a small imaginary component in the Young’s modu¬ 
lus of the material, with the magnitude adjusted for different 
Q m . The reflection spectra (Fig.[2}b) inset) show the expected 
Lorentzian dip as the waveguide mode frequency sweeps over 
the cavity, with the depth determined by the coupling rate rel¬ 


ative to the cavity intrinsic loss rate. For Q m = 1500 (typical 
for our devices), the reflection contrast is limited to ss 15 %, 
while an improved intrinsic Q m of 2.5 x 10 4 would enable crit¬ 
ical coupling. Alternately, the transition between waveguide 
and cavity, or the waveguide geometry itself, can be tailored 
to suitably increase the coupling rate. 

Experimentally, the IDT frequency response is measured 
by an electrical .S) 1 (reflection) measurement. The top panel 
of Fig. [2jd) shows the 5] 1 spectra for IDTs with pitch of the 
electrode fingers of 1.05 fim and 1.1 /rm. Optimal loading of 
the nanobeam’s localized mechanical mode by the RF drive 
occurs when the mode frequency lies within the IDT band¬ 
width. The bottom panel of Fig.[2]ri) shows the scaling of the 
IDT resonance frequency with finger pitch in simulation and 
experiment, along with the desired region of operation shaded 
in gray (corresponding to the nanobeam breathing mode fre¬ 
quency range). To ensure that some fabricated devices show 
the required spectral overlap, each IDT addresses an array of 
optomechanical cavities (Fig.[T]a)), and the nanobeam width 
is varied across the array, thereby tuning the localized me¬ 
chanical mode frequency with respect to a fixed IDT fre¬ 
quency 22 . 

While IDTs are a convenient means to generate propagat¬ 
ing acoustic waves for coupling to optomechanical cavities, 
the relatively small piezoelectric transduction coefficient of 
GaAs 2 ’ limits the efficiency with which acoustic energy can 
be converted to RF energy for electrical readout. Thus, for 
sensitive detection of the nanobeam optomechanical cavity’s 
motion when driven by the propagating acoustic wave, we use 
the strong photoelastic coupling for optical readout. The setup 
is shown in Fig.[2jc), where light is injected into the optome¬ 
chanical cavity using an optical fiber taper waveguide, and a 
vector network analyzer (VNA) drives the IDT with RF en¬ 
ergy and coherently detects the photocurrent signal transmit¬ 
ted past the optomechanical cavity, yielding an electro-optic 
521 measurement. The laser wavelength is positioned on the 
shoulder of the optical mode ( Q opt = 36600 ± 400, fit uncer- 
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FIG. 2: Coupling between propagating and localized phononic modes, (a) Finite element method (FEM) simulation of the surface acoustic 
wave generated by an input IDT (left) and its excitation of the propagating mode of a phononic crystal waveguide (right), (b) FEM simulation 
of the coupling between the phononic crystal waveguide and breathing mode of the nanobeam optomechanical crystal cavity. The inset shows 
the calculated acoustic reflection spectrum, for intrinsic mechanical quality factors of 1.5 xlO 3 (blue), 1.5x 10 4 (red), and 2.5 x 10 4 (green), 
respectively, (c) Experimental setup for optical readout of the nanobeam cavity coherent mechanical motion. The cavity is either driven 
through a phononic crystal waveguide excited by IDT 1, or through the phononic shield excited by IDT 2. APD = avalanche photodiode, (d) 
(Top) IDT reflection spectrum (Sn) for two different finger spacings. (Bottom) IDT resonance wavelength as a function of finger spacing, 
including both experimental results (blue circles) and FEM simulations (red squares), (e) (Top) Thermal noise spectrum of the nanobeam 
cavity’s localized mechanical mode. The y-axis span is 10 dB. (Bottom) Coherently detected motion of the nanobeam mechanical cavity as a 
function of RF frequency, both when the system is driven through a phononic crystal waveguide (red) and through the phononic shield (blue). 


tainty 95 % confidence intervals} 2 -^ so that phase fluctuations 
induced by motion are transduced into an intensity-modulated 
optical signal at the output. 

The bottom panel of Fig.|2|e) shows the measured electro¬ 
optic 521 in two instances. In the first case (red curve), a 
phononic crystal waveguide butt-coupled to the optomechani¬ 
cal cavity is adjacent to the RF-driven IDT (IDT 1, Fig.[2jc)), 
and a pronounced peak « 20 dB above the background is ob¬ 
served at the nanobeam cavity’s localized mechanical mode 
frequency, as confirmed by direct detection of its thermal 
noise spectrum (top panel of Fig. [2je)). In the second case 
(blue curve), an unperturbed phononic shield spans the re¬ 
gion between the cavity and RF-driven IDT (IDT 2, Fig. [2jc)). 
Here, the electro-optic 52i is suppressed by more than 30 dB 
relative to the first case, indicating that the phononic shield ef¬ 
fectively blocks the transmission of acoustic energy from the 
IDT to the cavity. 

Coherent Mechanical State Preparation and Readout 

Having demonstrated routing, we now further quantify the 
phonon injection process. Figure [3j a) shows the setup used, 
where a VNA excites a source IDT and reads out the coher¬ 
ent component of the photodetected light transmitted past the 
optomechanical cavity, with a phase shifter now added for 


phase control of the injected RF tone. In addition, a sec¬ 
ond RF signal generator drives an optical phase modulator for 
determining the optomechanical coupling rate go 18 -®, and a 
real-time spectrum analyzer (RSA) displays the total power 
spectral density (PSD), consisting of the incoherent thermal 
motion of the mechanical resonator, the coherent motion of 
the resonator due to coupling from the IDT-driven phononic 
waveguide, and the phase modulator calibration tone. 

Figure [3jb) displays the electro-optic S 21 measurement of 
coherently-driven motion similar to Fig. |3Je), but now for an 
optomechanical cavity connected to phononic waveguides at 
both its input and output (Fig. [Sj'a)). The results are similar, 
with the electro-optic 521 showing a peak corresponding to 
the cavity’s mechanical resonance frequency and which lies 
within the IDT’s resonance bandwidth. The importance of 
this overlap is shown in Fig. [3je)-(fj, which contrasts the RF 
power applied to the IDT to achieve a certain coherent intra¬ 
cavity phonon population when the localized nanobeam me¬ 
chanical mode is on-resonance (Fig. [3jc)) to when it is off- 
resonance with the IDT (Fig. [3jd)). For this comparison, 
two nanobeam optomechanical cavities from the same array 
were chosen, so that the IDT characteristics were fixed as the 
nanobeam mechanical frequency changed. Off resonance, the 
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FIG. 3: Acousto-optic modulation and coherent phonon detection, (a) Setup for determining the average coherent intracavity phonon 
number in the mechanical cavity. The coherent motion due to the RF drive is compared to the thermally-driven motion of the cavity. 
DUT = device under test, VNA = vector network analyzer, RSA = real-time spectrum analyzer, (b) IDT reflection spectrum (top) and optically- 
detected coherent mechanical motion (bottom) for an IDT and nanobeam that are spectrally aligned, (c)-(d) Thermal noise spectrum of the 
localized mechanical mode (red, RF turned off) and IDT electrical Si i spectrum (blue) for an IDT and nanobeam that are spectrally (c) aligned 
and (d) misaligned, (e)-(f) Corresponding total photodetected spectrum with the RF drive on. ApjJ 1 ’ = 43 ± 1.5 in (e). To achieve a similar 
= 50 ± 1.8 in (f), the RF drive power is increased by 20 dB. An additional tone due to phase modulation of the optical signal is used for 
calibrating the optomechanical coupling rate go- (g) Phase space diagram of the coherent mechanical motion for different values of the RF 
phase shifter. Each state has Ap£ h = 91 ± 3.3, and the VNA measurement bandwidth is 200 Hz. (h) (Left) Polar plot of (Vp(J h for different RF 

drive amplitudes. (Right) Total photodetected spectrum with the RF drive on and attenuated to produce (VpjJ h = 0.3 ±0.01. (i) Histogram of 
the the polar plot data from (h), shown as a function of motional amplitude, with ApjJ h and the measurement noise floor (obscured in the polar 
plot in (h)) indicated. Measurement bandwidth is 20 Hz. The uncertainty in the phonon number is due to the combined uncertainty in sample 
temperature, fit to the thermal noise spectrum, and RSA amplitude error and is a one standard deviation value 22 . The power spectral density 
plots in (e), (f), and (h) are referenced to a power of 1 mW (0 dB = 1 mW). 


RF power is increased by nearly 20 dB to achieve a similar 
strength of the coherent tone. 

We quantify the average coherent intracavity phonon num¬ 
ber by comparing the strength of the coherent tone with 
that due to the thermally-driven motion. That is: 

j\r c °h _ k^T S co h(Q coh ) 

ph s th (n,„) 

where k\>, is Boltzmann’s constant, T is temperature, £2„, is the 
mechanical mode frequency, 5 co h(^coh) is the power in the co¬ 
herent tone, and S t h(£ 2 ,„) is the integrated power in the thermal 
noise peak. A±° h = 43 ± 1.5 for Fig. ^e) and N^ h = 50 ± 1 .8 


for Fig. [3jf), quantifying the increased RF drive needed to 
achieve a given A±(J h when the nanobeam mechanical reso¬ 
nance and IDT are spectrally misaligned. Finally, we use the 
phase modulator to measure go/2n = 1.1 MHz ± 60 kHz, 
corresponding well with simulations®, and with a one stan¬ 
dard deviation uncertainty set by the uncertainty in the phase 
modulator V K 22 . 

Another way to characterize the device is to consider it as 
an acousto-optic phase modulator (the comparison is exact 
in the sideband-resolved regime, where the mechanical mode 
frequency exceeds the optical cavity linewidth) and determine 
its V n . For the device in Fig. [3jc), V n m 652 m\®. This can 
be reduced to V n < 100 mV by optimizing the path from the 
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IDT to the optomechanical cavity and minimizing reflections 
and scattering at interfaces. We note the qualitative difference 
between these acousto-optic modulators and previous demon¬ 
strations, in which the propagating acoustic waves were not 
coupled to localized mechanical resonances 1927 . Use of a lo¬ 
calized mechanical resonance significantly enhances the mod¬ 
ulation efficiency without creating a significant bandwidth re¬ 
striction beyond that already imposed by the IDT. 

To demonstrate coherent control of the resonator, we would 
like to transfer both the amplitude and phase of the applied RF 
voltage onto the cavity displacement. Figure[3jg) plots the in- 
phase and quadrature components of the photodetected signal 
on-resonance for a variety of different phases applied to the 
IDT. Here, the data is displayed as a function of average co¬ 
herent mechanical displacement, where the calibration proce¬ 
dure described in the previous paragraph is used to determine 
/Vpi" 1 ', and is converted to a displacement: 


n: 


coh 


ph 


1 

2 


f OCcav 

\ -^zpf 


2 


( 2 ) 


where a ca v is the magnitude of the cavity displacement 
and x zp f = \Jn m i s th e zero-point motional amplitude 
(m e ff ~ 0.5 pg is the simulated effective mass). The data 
shows that, by varying the amplitude and phase of the ap¬ 
plied RF signal, we can imprint a corresponding amplitude 
and phase on the localized mechanical mode, and thus initial¬ 
ize it at an arbitrary location in phase space. These measure¬ 
ments confirm that coherence is preserved through the process 
of converting an applied RF voltage to an acoustic wave that 
propagates through a phononic crystal waveguide and couples 
to the localized mode of the optomechanical crystal cavity. By 
injecting propagating acoustic waves with similar amplitudes 
and varying phase difference through two phononic waveg¬ 
uides that excite the same optomechanical cavity, we can opti¬ 
cally read out constructive and destructive interference effects 
occurring within the mechanical resonator 22 . 

The combination of coherent detection and the large op¬ 
tomechanical coupling strength enables measurement of a 
weak coherent intracavity phonon population on top of the 
~ 2600 thermal phonons in the cavity. In Fig. [3j h), we re¬ 
duce the amplitude of the RF drive and measure the mechan¬ 
ical motion on-resonance, displayed as a polar plot in terms 
of the phonon number and phase angle. Figure |3ji) shows the 
data plotted against motional amplitude for the three smallest 
RF drive amplitudes, where for each the data has been his- 
togrammed into 15 equally spaced bins. /\'j|j >h « 1.1 is well- 
separated from the measurement noise floor (RF drive set to 
zero), and as low as « 0.1 is resolvable. We also display 
the PSD for a measurement with A^j 11 < 1 (Fig.plh) inset). 
Coherent Interactions between RF, Optics, ana Mechanics 
The focus of our work so far has been establishing phononic 
crystal waveguides that are sourced by IDTs and effectively 
couple acoustic energy to the localized mechanical mode of a 
nanobeam optomechanical cavity. In these experiments, large 
go enables sensitive optical readout of the mechanical motion, 
but was not used to manipulate the mechanical resonator. We 


next consider experiments in which the mechanical resonator 
is simultaneously excited using both RF signals (via IDTs and 
phononic crystal waveguides) and optical fields. This high¬ 
lights a key feature of the GaAs optomechanical circuit plat¬ 
form, which is that the motion of the mechanical resonator can 
be manipulated with similar facility along either the optical or 
RF channel. 

Figure |4|a) shows a schematic level diagram indicating the 
couplings and decay channels of the different components in 
our system. If we ignore the states on the right hand side of 
the diagram that depict the phononic waveguide, its feeding 
from an IDT, and coupling to the phononic cavity, we recover 
the characteristic A-system configuration in which optome- 
chanically induced transparency (OMIT; 17 S 28 has been ob¬ 
served in a wide variety of systems. Adding the new chan¬ 
nels provides us with another knob to coherently control and 
probe the system properties and observe new phenomena. The 
interference effect we describe in this section occurs in the 
phononic cavity, which can be driven optically via the beat¬ 
ing of the probe and control lasers that are detuned from each 
other by the localized mechanical mode frequency (as in typ¬ 
ical OMIT), or acoustically through the phononic waveguide. 
These two pathways can be tuned to have equal amplitudes 
and opposite phases, resulting in a cancellation of the coherent 
motion of the mechanical resonator. This effect is analogous 
to coherent population trapping in atomic systems 29 30 . 

Starting with the coupled equations of cavity optomechan- 
ics, one can derive an expression for the cavity displacement 
amplitude p,/ + (see supplementary material 2 ): 


(W 


-ig 0 (a 0 a* d _) - 

i{&m ~ &d) + \ 


(3) 


The first term in the numerator corresponds to the optical 
drive term due to beating of the intracavity control (ao) and 
probe (a,/ ) beams and the second term corresponds to the 
acoustic drive term (P;«, o^) due to acoustic waves coupling 
to the cavity (intrinsic decay rate y ’,) through the phononic 
waveguide (coupling rate y e ). Appropriate choice of the am¬ 
plitude (Pi„ o) and phase (\|/) of the acoustic drive term pro¬ 
duces the condition for cancellation of the mechanical mode 

(Prf+ = 0): 



P,„.oe' V = 


TSoaoaJ, 


— i(A + Q.d) + y 


(4) 


wherein 0C, nr represents the input probe signal incident on 
the cavity, A the laser-cavity detuning, £2 ( / the probe beam de¬ 
tuning from the control beam frequency, K e the optical cavity 
coupling rate, and K, the intrinsic optical cavity decay rate. 

When this condition is satisfied, the coherent cavity dis¬ 
placement tends to zero and the cavity transmission is fiat 
over the mechanical resonator’s bandwidth (the optomechan¬ 
ical cavity acts as a pure optical cavity to the probe signal). 
The mechanical cavity effectively enters a “dark” state for 
phonons analogous to the dark state for photons in coher¬ 
ent population trapping. It is important to note that only the 
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FIG. 4: Acoustic wave interference analogous to coherent population trapping, (a) Schematic level diagram indicating how the localized 
phononic cavity is populated either through the optical channel (via photoelastic coupling to the photonic cavity) or through RF-generated 
propagating acoustic waves. These two pathways can coherently interfere and cancel the mechanical motion, in an effect analogous to atomic 
coherent population trapping, (b) Experimental setup. One part of the VNA output is sent to an optical phase modulator, after which interfer¬ 
ence between the generated optical sideband and carrier drives the phononic cavity optically. The other part of the VNA output is amplified 
and sent through an RF phase shifter into the input IDT to generate a propagating acoustic wave to couple to the phononic cavity. Readout 
of the cavity’s coherent mechanical motion is done optically, (c) Experimental data (left) and theoretical calculations (center) of the coherent 
mechanical motion as a function of RF drive frequency and for increasing RF power. Pure optomechanically-driven motion (RF = 0) is shown 
in the black curves, acoustic coherent population trapping is shown in cyan, and RF-dominated motion is shown in dark gray. The rightmost 
set of curves shows the corresponding theoretical calculations for the cavity displacement amplitude, (d) Phase space diagram for different 
values of the RF amplitude and phase. The black dashed semicircle is the far off-resonance background. The four numbered points, with 
frequency response curves shown in (e), correspond to (i) pure optomechanically-driven motion, (ii) acoustic coherent population trapping, 
(iii) RF-dominated motion, and (iv) Fano lineshape behavior. The measurement bandwidth in (d)-(e) is 200 Hz. 


coherently driven component of the cavity displacement be¬ 
comes zero, as the cavity is still driven by incoherent thermal 
motion and has a finite power spectral density. We also em¬ 
phasize the distinction between this effect and other interfer¬ 
ence phenomena observed in cavity optomechanical systems. 
OMIT is an optical interference effect in which probe pho¬ 
tons interfere with photons scattered by the driven mechanical 
resonator. A closely-related effect in piezo-optomechanical 
systems has been termed electromechanically-induced optical 
transparency 13 or microwave-assisted transparency 14 , where 
interference still occurs in the optical domain, but the me¬ 
chanical resonator is driven by an RF voltage. In these sit¬ 
uations, which we also observe in our system, the mechanical 
resonator has a finite coherent displacement amplitude at all 
times. In contrast, during acoustic coherent population trap¬ 
ping, the mechanical resonator’s coherent displacement am¬ 
plitude is zero. Coherent wavelength conversion between two 
optical modes coupled to the same mechanical mode 31 33 is 
similar in that it can also be viewed as being mediated by an 
optomechanical dark mode 32 ®, although coherent mechani¬ 
cal motion is suppressed purely though optomechanical cou¬ 


pling. 

Figure |4jb) shows the experimental setup used to probe the 
coherent interaction between RF, optics, and mechanics in our 
system. The RF output from port 1 of a VNA is split in two 
and drives both an electro-optic phase modulator and the IDT. 
An RF amplifier and phase shifter control the amplitude and 
phase of the RF signal driving the IDT with respect to that 
driving the phase modulator, and photodetected light trans¬ 
mitted past the optomechanical cavity is detected by port 2 of 
the VNA. The control laser power is fixed and is blue-detuned 
by a mechanical frequency from the optical cavity resonance. 

Figure |4]c) shows a series of electro-optical .SAi measure¬ 
ments performed as the VNA output frequency is swept across 
the mechanical resonance frequency and the RF power into 
the IDT is increased, along with calculations of the expected 
S 21 spectrum and corresponding mechanical resonator dis¬ 
placement amplitude 22 . The topmost black curves are for the 
IDT channel turned off, and correspond to the pure OMIT 
case. We then increase the RF power while keeping the 
phase fixed to satisfy the condition for acoustic wave inter¬ 
ference. As the RF power is increased, the OMIT peak is re- 
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duced and at a particular RF power (cyan curves in Fig. |4jd)), 
the transmission spectrum reaches zero on-resonance and is 
approximately (perfectly) spectrally flat in experiment (the¬ 
ory). This corresponds to complete suppression of the co¬ 
herent cavity displacement and removal of the transparency 
window that was initially induced by the optomechanical in¬ 
teraction. The inability to achieve perfect spectral flatness in 
experiment is attributed to the actuation of other mechani¬ 
cal modes that occurs electrically but not optically 22 . If we 
further increase the RF power, the cavity is primarily driven 
by the phononic waveguide channel and the cavity displace¬ 
ment amplitude increases correspondingly (magenta, yellow, 
dark gray curves). Since the phase is set for destructive in¬ 
terference, the probe transmission drops and for these higher 
RF powers, the interference feature within the optical cav¬ 
ity spectrum re-appears (the system is now in the regime of 
electromechanically-induced optical transparency 13 14 ). Here, 
the mechanical breathing mode is in the same state as in the 
pure OMIT case, except that the displacement has the oppo¬ 
site phase and so transparency is converted to absorption. 

A more complete picture is observed by measuring the in- 
phase and quadrature components of the on-resonance pho- 
todetected signal (Fig. |4jd)). Starting with the pure OMIT 
condition (IDT turned off. Fig. |4|e), curve i), as we increase 
the RF drive to the IDT, for the right choice of phase the cavity 
displacement approaches zero (Fig.|4je), curve ii), and the sig¬ 
nal on-resonance is equal to its off-resonance value (Fig. |4je), 
dashed circle). Further increase in the RF power moves the 
system into the region of induced absorption (Fig. |4je), curve 
iii), while other parts of phase space are accessed by vary¬ 


ing the phase applied to the IDT. For an arbitrary phase, the 
response of the system is given by an asymmetric Fano line- 
shape (for example, Fig.^e), curve iv). 

Conclusions 

We have presented an integrated piezo-optomechanical cir¬ 
cuits platform in which fiber-coupled GaAs optomechanical 
crystal cavities are interfaced with phononic waveguides and 
interdigitated transducers. This system combines localized 
and interacting 1550 nm photons and 2.4 GHz phonons with 
excitation and readout along either the optical or mechanical 
channel, with the latter connected to the RF domain via IDTs. 
A distinguishing feature of GaAs is that it is both piezoelec¬ 
tric and has a strong photoelastic effect, so that the mechanical 
resonator can be effectively manipulated with equal facility 
by either RF or optical fields. Going forward, these capabil¬ 
ities may enable the storage and retrieval of either optical or 
acoustic pulses, while phononic routing between cavity op¬ 
tomechanical nodes has been envisioned for both classical and 
quantum device applications 35 36 . Further technical develop¬ 
ment may include optimized IDTs 23 37 for increased acous¬ 
tic wave generation efficiency, or use of stimulated Brillouin 
scattering as an optical mechanism for acoustic wave genera- 
tioiPS Developments in chip-based phononics 38 39 39 42 can 
inform future phononic circuitry to interface with the optome¬ 
chanical cavities, while the incorporation of InAs/GaAs quan¬ 
tum dotPES 

which have been used in a number of quantum 
optics application^ anc j are fully compatible with the device 
fabrication described here, would add qualitatively new func¬ 
tionality. 
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SUPPLEMENTARY INFORMATION 

I. FABRICATION PROCEDURE 


The epitaxial material used in this work consists of a 220 nm thick GaAs layer on a 1.5 /./m thick Alo. 7 Gao. 3 As sacrificial 
layer. The samples were spin-coated with positive tone electron beam resist and baked at 180°C for 120 s. The IDT patterns 
were exposed in a 100 keV direct write electron beam lithography system with a beam current of 2 nA and nominal dose of 
580 ,uC/cnr. After exposure, the electron beam resist was developed using MIBK:IPA 1:3 solution (90 s). To remove residual 
undeveloped resist, the samples were exposed to an O 2 plasma (6.7 Pa = 50 mTorr, 75 W) for 5 s. The IDT metallization was 
carried out in an electron beam evaporator with the metal stack Cr (5 nm)/Pt (15 nm)/Au (30 nm) deposited in succession. The 
lift-off was carried out by soaking the samples overnight in acetone and using gentle sonication. 

For the nanobeam overlay, the samples were spin-coated with positive tone electron beam resist and baked at 180°C for 120 s. 
The nanobeam patterns were exposed in a 100 keV direct write electron beam lithography system with a beam current of 200 
pA and nominal dose of 250 ,uC/cm 2 . The nanobeam patterns were aligned to the IDT using alignment marks that were exposed 
along with the IDTs. After exposure, the electron beam resist was developed using hexyl acetate (65 s). To remove residual 
undeveloped resist, the samples were exposed to an O 2 plasma (6.7 Pa = 50 mTorr, 75 W) for 10 s. The nanobeam patterns were 
then transferred to the underlying GaAs layer using an inductively coupled plasma reactive ion etcher with an Ar/CL chemistry. 
The electron beam resist was stripped using trichloroethylene, and the nanobeams were undercut with a timed wet etch using 
either concentrated (49 %) HF solution or (NFLO 2 S and dilute HF (50:1 volume dilution in FLO). 


II. SAW RESONANCE SIMULATION 


To calculate the resonance frequency for surface acoustic waves for our material stack, we model a one wavelength width 
cross-section of the stack with periodic boundary conditions and solve for the resonance frequency of the stack. We used 
the following parameters to model Al v Gai^ A As with x = 0.7: = — 0.16 — 0.065x C/m 2 , Cn = (118.8+ 1 Ax) GPa, 

C 12 = (53.8 + 3.2*) GPa and C 44 = (59.4 — 0.5x) GP^. The calculated frequencies are uniformly lower than the measured 
frequencies, which we attribute to imperfect knowledge of the density and thickness of the metallic stack (Cr/Pt/Au) comprising 
the IDT. 


III. FREQUENCY RESPONSE ANALYSIS FOR ACOUSTIC MODE COUPLING 

For the simulation results shown in Figure 2(b) in the main text, the surface acoustic wave was launched by applying a 
prescribed displacement of 1 nm (the simulation was linear in displacement, the starting displacement amplitude is arbitrary) 
in the x-direction to the end facet of the line-defect waveguide (x is the direction of wave propagation). The x-displacement is 
chosen so as to launch a symmetric Lamb wave in the defect waveguide and nanobeam that has the right ^-symmetry to couple 
to the nanobeam breathing mode. 

To calculate the coupling efficiency shown in Figure 2(b) we define the acoustic Poynting vector (Pj) 

P, = ] 2 Re(-T u vr) (SI) 

where 7)y is the stress and v, is the particle velocity. The reflection spectrum is calculated by doing two simulations, one with the 
nanobeam cavity and a second with a bare waveguide whose width is equal to that of the nanobeam cavity. The Poynting vector 
is calculated at the input port of the cavity and the reflected power is defined as: 


. Grans.cav 

K = 1 - 

Grans,wvg 

where the transmitted powers (Grans.cav and Grans,wvg) are all calculated in the direction of the propagating mode ( x). 


(S2) 


IV. NANOBEAM BREATHING MODE FREQUENCY DEPENDENCE ON WIDTH 

To ensure that the resonance frequency of the nanobeam optomechanical cavity’s mechanical breathing mode lies within the 
excitation bandwidth of the IDT, we excite an array of cavities with varying beam width with the same IDT. Figure [87]+) shows 
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the dependence of the breathing mode frequency on the width of the nanobeam and indicates that a 40 nm tuning in the width 
results in an approximately 155 MHz tuning in breathing mode frequency. In fabricated devices, the width of the individual 
beams within the nanobeam array is varied in 5 nm steps. 




FIG. SI: (a) Calculated nanobeam optomechanical crystal cavity mechanical breathing mode frequency as a function of beam width, (b) 
Schematic of the unit cell of the phononic shield showing the different dimensions. 


V. PHONONIC SHIELD BAND STRUCTURE 


The dispersion diagram for acoustic wave propagation in the T-X direction for the phononic shield (whose unit cell is shown 
in Figure [ST[ b)) is shown in Figure [S2] The modes of a unit cell with periodic boundary conditions in the propagation direction 
(r-X) were calculated using a finite-element solver to determine the dispersion diagram. The modes are classified as odd / even 
according to the symmetry in the out of plane (z) direction. GaAs is modelled as an orthotropic elastic material with parameters 
E x = 121.2 GPa, E y = 121.2 GPa, E z = 85.9 GPa, V lv = 0.0209, \ yz = 0.4434, M xz = 0.312, G xy = 32.5 GPa, G yz = 59.4 GPa, and 
G xz = 59.4 GPa. 


(a) Even (b) Odd 



FIG. S2: Calculated phononic bandstructure for acoustic waves with (a) even and (b) odd symmetry in the out-of-plane direction (z), propa¬ 
gating along the T-X direction in the phononic shield. The presence of a complete bandgap for the frequencies of interest (2.25 GHz to 2.75 
GHz) is clearly seen. 


We note that a wide variety of phononic shield structures based on circular / square holes have been studied 2 5 . For the GaAs 
material system, achieving a complete bandgap within such geometries necessitates high hole filling fractions, and in the 2.4 
GHz frequency range this requires patterning dimensions at the sub-25 nm length scale, making fabrication challenging. The 
advantage of the cross structure is that it allows one to achieve a complete acoustic band-gap with minimum feature size of 100 
nm, which can be easily obtained by standard electron-beam lithography. 
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VI. BASIC OPTOMECHANICAL SPECTROSCOPY 

Figure |S3[ a) shows the optical transmission spectrum of the device in Figure 3(g)-(h) and Figure 4 (main text) whose 
Qopt = 36600 ± 400 (uncertainty from the 95 % confidence interval of a Lorentzian fit to the data). Figure [S3f b) shows 
the corresponding thermal noise spectrum for the localized mechanical mode with Q m = 1400 ± 5 (uncertainty from the 95 % 
confidence interval of a Lorentzian fit to the data). The optomechanical coupling rate go/2 k for the device was measured to be 
1.1 MHz ± 60 kHz using the procedure discussed below. 



FIG. S3: (a) Transmission spectrum of the cavity optical mode, (b) Thermal noise power spectral density (PSD) of the nanobeam breathing 
mode (as measured on the RSA). Blue curves represent experimental data, Red: Lorentzian fits. The uncertainty in the optical and mechanical 
quality factors represent Lorentzian fit uncertainty, 95 % confidence intervals. 


VII. g 0 CALIBRATION 

We use a phase modulator to calibrate the optomechanical coupling rate The basic idea is to relate the modulation pro¬ 
duced by sending light through the cavity optomechanical system (which is driven by its contact with the thermal environment) 
to a direct phase modulation applied with an electro-optic phase modulator. Because both undergo the same transduction func¬ 
tion, the ratio of the integrated power in their photocurrent RF spectra will be related to the ratio of go and the phase modulator’s 
modulation index. We go through this derivation below. 

The root-mean-square (rms) thermal displacement amplitude (a t h e rmai) is related to the temperature (T) using the equipartition 
theorem: 


1 , , 1 

2 m efA"n “thermal = ^ T ( S 3 ) 

where m e ff is the motional mass of the resonator and is the mechanical mode frequency. We define a thermal modulation 
index: 


P thermal — 


^thermal&om 


where the optomechanical coupling parameter (gom) is defined by: 

co(a) = ©0+gomCt 

with CDo being the unperturbed cavity frequency. This leads to: 


o2 _ T <g om 


(S4) 


(S5) 


(S6) 


The optomechanical coupling rate go is: 
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go — gom-^zpf (S7) 

where the (rms) amplitude of the zero point fluctuation is: 


Xzpf -/ 2 mln m 


(S8) 


which can be directly calculated as y/ (0|x 2 |0) from the ground state wavefunction |0) of the simple harmonic resonator. 
We then relate go to pthcrmaf 


Pthermal — ^3 #0 


The modulation index (p pm ) of the phase modulator is defined as: 


(S9) 


Ppm = (S10) 

me 

where y s j g is the signal amplitude and V K is the modulator half wave voltage. 

By comparing the (integrated) powers in the cavity mechanical mode signal (S cav (f2 m )) and the phase modulator signal 
(S'pmfflpm)) obtained from the electronic spectrum analyzer, we get: 


2 _ fifim ^2 o2 ^cav(^m) 
t?0 “ 2k B T mP P m S pm (£ 2 pm ) 


(Sll) 


VIII. EXTRACTING CAVITY DISPLACEMENT AND PHONON NUMBER DUE TO EXCITATION VIA THE PHONONIC 

WAVEGUIDE 

The modulation index due to the propagating wave p saw is given by: 


Psaw 


O^cavgom 


Q 


saw 


(S12) 


wherein a ca v is the cavity displacement (rms) due to the coupling from the propagating acoustic wave. 

To extract the average coherent intracavity phonon number A^^ h , we compare the area under the coherent acoustic wave peak 
with the integrated area of the mechanical mode thermal noise spectrum, noting that the average thermal phonon occupation 
number is given by N t /, ss k B T / 


Af c °h _ k B T Scoh(Qcoh) 

ph s th (kl m ) 


(S13) 


We can now relate the cavity’s coherent rms displacement amplitude a ca v to N c p ^ using: 


N^ h hn m = m et [Q.l a 2 cay 


which can be re-written as: 


(S14) 


N c S h = - 

ph 9 


Ctcav \ 
•*zpf / 


(S15) 


where we have substituted x zp f = 


h 

2m e ff£l m ' 
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IX. ACOUSTO-OPTIC MODULATOR V K CALCULATION 


To determine the equivalent V n of our piezo-optomechanical circuits when considering them to operate as an acousto-optic 
phase modulator, we use the relationship: 


T j ksaw 

^ = — fS 16) 

Psaw 

To determine V saw , the applied RF voltage that is converted to a propagating surface acoustic wave, we use the resonant dip in 
the Sn spectrum of our IDT to determine the transmitted RF power (in our case -0.25 dB) and assume a 50 Q load to determine 
the corresponding voltage. The modulation index due to the SAW is extracted using the procedure described in the previous 
section. 


X. ACOUSTIC WAVE INTERFERENCE FROM TWO IDTS 


The main text discusses acoustic wave interference that occurs when the mechanical mode is driven optically and electrically 
with equal amplitude and opposite phase. The same phenomenon occurs when the mechanical resonator is pumped by two IDTs 
through two phononic waveguides with equal amplitudes and opposite phase. To demonstrate this, we use the device shown in 
Figure 3(a) in the main text, and modify the experimental setup of Figure 2(c) so that the output of the VNA is split in two, 
with each path fed to its own RF amplifier and then connected to an IDT. The resulting data is presented in Figure [S4] The two 
IDTs when turned on alone (ii and iii in the polar plot) result in a regular electro-optic S 21 measurement as in Figure 2(e), for 
example. The asymmetry between the spectral responses can be attributed to fabrication-induced asymmetry in the device and 
the asymmetric coupling induced due to the fiber taper touching down on the device. When the two IDTs are both fed by an 
RF signal with a varying phase difference, one can sweep out a circular trajectory in phase space. The intersection of this circle 
with the origin corresponds to the acoustic wave cancellation condition. The mechanical resonator is loaded in a dark state by 
the destructive interference between the excitation amplitudes of the two IDTs. For an arbitrary phase difference between the 
two IDTs, the response in general has an asymmetric Fano-like lineshape. 




FIG. S4: Phase space diagram for different values of the RF amplitude and phase for a nanobeam optomechanical crystal cavity pumped by 
two phononic waveguides, each of which is sourced by its own IDT. The six numbered points, with frequency response curves shown in the 
side panel, correspond to (i) acoustic wave cancellation, (ii) IDT 1 ON, IDT 2 OFF, (iii) IDT 1 OFF, IDT 2 ON, and (iv, v and vi) Fano 
lineshape behavior for varying phase difference between IDT 1 and IDT 2. The measurement bandwidth in these experiments is 200 Hz. 
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XI. ACOUSTIC COHERENT POPULATION TRAPPING THEORY 

The equations of motion for the intracavity optical field amplitude and the cavity mechanical displacement amplitude are: 


a = — (iA + -y-)a — igoa(b + b^) 



CS17) 


b=-(iQ. m +^)b-ig 0 a f a -y (S18) 

with a(V) representing the annihilation (creation) operator for the intracavity optical field, bib') the annihilation (creation) 
operator for the mechanical displacement. A represents the detuning of the control beam from the optical cavity resonance 
frequency, K, the intrinsic decay rate of the optical cavity, and K e is the extrinsic decay rate (coupling rate) to the waveguide. 

represents the mechanical mode frequency, y, the intrinsic decay rate of the mechanical cavity and y e is the extrinsic decay 
rate (coupling rate) to the phononic waveguide. a lN and /?,„ represent the optical and acoustic field strengths in the photonic and 
phononic waveguides respectively, go represents the vacuum optomechanical coupling rate. 

We make the following ansatz: 


a => a = ao + v.d+e + a d~e l£ldt 


(S19) 


P = Po + P«/+«““ V (S20) 

Physically, this ansatz correponds to linearizing around the steady state photon (ao) and phonon (po) amplitudes. For the 
mechanical motion: 


IV- = p; /+ 

Substituting equations S19 and S20 into equation [ST8] and collecting terms at DC: 


(S21) 


p 0 = zv°ni 2 




(S22) 


The true numerator is —/go(|dCo| 2 + |oc^/+| 2 + |a ( /_| 2 ) but we neglect the latter two terms under the assumption that |ao| 2 >> 

| a </+| 2 + | ok /-| 2 . 

We can write the incident acoustic displacement as: 


P,„ = P i, 


fie 


'¥ 


(S23) 


wherein p„, o is the amplitude of the incident acoustic wave and \|/ represents the (variable) phase difference between the RF 
signal applied to the p hase mo dulato r (PM) and t he ID T. 

and collecting terms at e~ ,Q,dt : 


Substituting equations 


S19 


and 


S20 


into equation 


S18 


IV+ 


-ig 0 (a 0 a* d _+a d+ a* 0 ) - p,„,oe ,v 

i(&m ~ Ql) + \ 


(S24) 


We can now see that the mechanical resonator’s coherent amplitude p,/ | is driven optomechanically, through beating of the pump 
(ao) and probe sidebands (a d + an d a,/ ), as well as through the phononic channel by p m . When the amplitude of these two 
terms is equal, the phase t|/ can be tuned to achieve cancellation with p^ + =0, corresponding to the acoustic coherent population 
trapping condition. 

Having determined an equation for the mechanical resonator’s coherent displacement amplitude, we now consider the system’s 
optical response. The input field can be written as: 
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OCm = CLinfl (1 + ^^ + i^e ^) 


(S25) 


where (|) is the modulation index that represents the action of the optical phase modulator on the input optical pump field. 
Substituting equations |S19|and|S20|into equation [STT] and collecting terms at DC: 


a 0 = -t 


Y a in, 0 


z'(A + 2go(Po + Po))' 


Substituting equations S19 and S20 into equation S17 and collecting terms at e 




(S26) 


Ud+ = 


2igoOCoP</+ Y 2 ^ VUL 
i( A — fid) + y 


Substituting equations S19 and S20 into equation S17 and collecting terms at e 




(S27) 


- 2 zgoOC 0 pj + ~ \/^a, 


ad- = 


i(A + Q. d ) + y 


The two sidebands at the output are: 


(S28) 


ttd - .oul — T \l 2 ^d+| 


(S29) 


\a d ,oui\ — | a ddl 


-a d - 


(S30) 


In the resolved sideband case (when the pump is blue-detuned), we can derive a condition for observing acoustic wave 
interference. The relevant equations, in the limit that only the a,/ sideband and p f ; ( matter, are: 


Pd+ = 


-igo(a 0 a* d _) - Y/fp^oe'V 
*(fi m — fid) + T 


(S31) 


2z#0«0p rf _ Y 2 ain,— 
i (A + Q . d ) + y 

Acoustic coherent population trapping corresponds to P f / + = 0. In this limit, a. d reduces to: 


(S32) 


a d - 



i (A + H d ) + y 


(S33) 


If we substitute this back in equation |S31| we get the condition for the amplitude and phase of the propagating acoustic wave 
input to engineer acoustic coherent population trapping: 



Pm ! o<?' v = —igoaoa* d _ 


Substituting a d - from equation S32 leads us to: 



-iyftgoaoa’m,-- 

— z(A + fi f /) + y 


(S34) 


(S35) 
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XII. COMPARISON BETWEEN OMIT, EMIT, AND ACOUSTIC CPT 


A number of different coherent interference effects have been observed in cavity optomechanical systems, and more recently 
j n piezoelectric ally actuated cavity optomechanical systems 11 12 . Here, we compare these different effects with those ob¬ 
served in our experiments. A key conclusion is that while these previous demonstrations have involved interference in the optical 
domain, the acoustic coherent population trapping effect we have observed is a novel cancellation occuring in the mechanical 
domain. 

Figure [S5] schematic ally depicts four different physical situations: (a) waveguide-cavity coupling in a system without any op¬ 
tomechanical interaction; (b) optomechanically-induced transparency 8 "^, (c) electromechanically-induced optical transparency 
(EMIT) or microwave-assisted OMIT 11 12 ; and (d) acoustic coherent population trapping. The pump is detuned one mechan¬ 
ical resonance away from the optical cavity line-center and we monitor the probe transmission by looking at the coherent RF 
spectrum of the photodetected signal (generated by the beating between the pump and the probe) transmitted past the cavity. We 
assume that the system is sideband-resolved, so that only one of the sidebands produced by the phase modulator interacts with 
the cavity. Our system is sufficiently close to the sideband-resolved regime (K„ p; /27t « 5.2 GHz, G m /2n ss 2.4 GHz) for this 
assumption to be largely valid. In the calculations presented in the main text, the influence of both sidebands was included to 
model the system as closely as possible. 

Starting with a simple side-coupled waveguide cavity system (assuming no optomechanical coupling), an incident optical 
probe will show a characteristic Lorentzian dip in its transmission spectrum corresponding to the optical cavity resonance 
(shown in Figure [S5f a)). If we now turn the optomechanical coupling on, when the frequency difference between the pump 
and the probe approaches the mechanical resonance frequency, their interference drives the mechanical motion of the cavity and 
scatters photons from the pump to the probe beam frequency. These additional probe photons can interfere either constructively 
or destructively (depending on the pump detuning) with the transmitted probe beam (from the regular waveguide-cavity scenario) 
and leads to the opening of a transparency window in the transmission spectrum (Figure [S5| 'h)), a process commonly referred to 
as optomechanically induced transparency (OMIT) 9 !0 . It is important to keep in mind that the interference occurs in the optical 
domain and the mechanical mode is coherently driven (has a non-zero displacement amplitude). 

To observe OMIT, the optomechanical coupling rate go/ln has to be sufficiently strong so that the beating between the pump 
and probe can drive the mechanical motion with sufficient amplitude to scatter photons from the pump to the probe beam 
frequency. In systems with small go/2tt, the mechanical motion can be coherently driven by a different mechanism (usually by 
acoustic waves using the piezoelectric effect) and as long as the RF signal driving the mechanics is derived from the same source 
as the signal driving the electro-optic phase modulator, the photons scattered by the resonator from the pump will be phase 

c)). This has been referred to as 

mm 


coherent with the probe photons and one can see the induced transparency effect (Figure S5 
electromechanically induced transparency or microwave-assisted transparency in the literature 


In effect, it is not different 

from OMIT, except the source driving the coherent mechanical motion is electrical rather than optical, and the character of the 
interference (constructive, destructive, or Fano-like) is set by the relative phase between the electrical drive and optical probe. 

Finally, in a system with both strong go/2tt and electrical drive, as we have presented in this work, the mechanical resonator 
can be driven both optically (due to the beating between pump and probe, like in OMIT) or electrically (like in EMIT). By 
choosing the amplitude and the phase of the electrical drive, one can see either constructive or destructive interference in the 
mechanical domain (Figure |S5f d)). Unlike pure OMIT / EMIT, the interference actually occurs in the acoustic domain and 
in case of destructive interference, the coherent component of the mechanical motion is zero. In this case, the system shows 
coherent population trapping for phonons analogous to the CPT observed in atomic systems. One can also say that the system 
has reached a mechanically ’dark’ state, and the transmission goes back to that of a bare optical waveguide-cavity system. 

The term ’coherent population trapping’ was chosen because the overall system, schematically depicted in Figure 4a and 
re-displayed in Figure |S6[ a), can be viewed from the perspective of the phononic cavity as an effective three-level system, in 
which the phononic cavity is populated through either an optomechanically-mediated path or an RF-driven path (Figure |S6f b)). 
The mechanical dark state - cancellation of the phononic cavity’s coherent mechanical motion - occurs when the strength of the 
two transitions is equal (and phase is opposite). This is analogous to the situation in atomic physics, where coherent population 
trapping refers to an interference occuring when the two terms are equal in amplitude, while EIT is used when the control field 
is much stronger than the probe field *2 .ST 

Finally, we note that an effect analogous to OMIT, but now consisting of a transparency in the spectrum of propagating 
acoustic probe phonons due to coupling to intracavity photons, should also be observable in this system. 


XIII. ABSENCE OF PERFECT CANCELLATION IN COHERENT POPULATION TRAPPING 


The inability to achieve perfect cancellation (a flat transmission spectrum in either the cyan curve from Figure 4(c), or Figure 
4(e), curve ii) can be attributed to the difference between the optomechanical (K om ) and electromechanical (K sm ) transduction 
spectra for different mechanical modes in the optomechanical cavity. To illustrate this, we consider the difference between 
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FIG. S5: Illustration of different interference effects observable in the piezo-optomechanical devices, (a) Standard optical waveguide- 
cavity coupling, in which the mechanical mode does not play a role (optomechanical and electromechanical coupling are set to zero); (b) 
optomechanically-induced transparency (OMIT); (c) electromechanically-induced optical transparency (EMIT); (d) acoustic coherent pop¬ 
ulation trapping (CPT). For simplicity, the figures consider the sideband-resolved regime in which only one sideband produced by phase 
modulation of the optical pump interacts with the system. 
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FIG. S6: (a) Schematic level diagram for the piezo-optomechanical circuit, re-displayed from Figure 4(a) of the main text, (b) Effective 3-level 
A system for phonons indicating the optomechanically-driven and RF-driven pathways. 


the OMIT response wherein we are primarily measuring the optomechanical transduction coefficient and the electro-optic Si 1 
measurement wherein we measure the product of the two: 


PSD om ,t - K om (Cl) (S36) 


PSDeo-sii « K em (Fl)K om (Fl) (S37) 

The optomechanical transduction (K om ) is calculated by an overlap integral between the localized optical and mechanical 
modes whereas the electromechanical transduction is estimated by an overlap integral between the localized mechanical mode 
and the surface acoustic wave incident on the cavity. When a surface acoustic wave is incident on the cavity, it excites a 
superposition of all the mechanical modes of the cavity which lie within the given frequency range and the corresponding mode 
amplitudes are given by the corresponding transduction coefficients. 



FIG. S7: (a) Normalized probe transmission for a nanobeam optomechanical cavity when the mechanical mode is coherently driven (i) optically 
(OMIT, green curve) and (ii) acoustically, through a surface acoustic wave excited by the IDT (SAW, blue). 

This helps us understand why the two spectra appear different (shown in Figure [S7] (a)). In an OMIT measurement, the 
only mode that is excited with significant amplitude is a mode with strong optomechanical coupling (in this case, the localized 
breathing mode) and hence, the power spectral density has a Lorentzian lineshape corresponding to a single mode. When a 
propagating acoustic wave is incident on the cavity, other modes which might have small optomechanical coupling but sig¬ 
nificant electromechanical coupling (due to a strong overlap with the propagating acoustic wave) also appear in the measured 
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transmission spectrum. Since the acoustic wave interference condition requires both the modes to have similar optomechanical 
and electromechanical transduction for the cancellation to occur, only the breathing mode is cancelled while the other excita¬ 
tions are not. This leads to the appearance of residual peaks and a non-flat background in the transmission spectrum. (Note: we 
believe these additional mechanical modes occur due to fabrication imperfections in the dimensions of the elliptical holes in the 
nanobeam cavity). 


XIV. ACOUSTIC WAVE INTERFERENCE FOR A NANOBEAM CAVITY WITH TWO CLOSELY SPACED MECHANICAL 

MODES 


As illustrated in the previous section, different closely spaced mechanical modes can have vastly different optomechanical and 
electromechanical coupling coefficients. This enables us to demonstrate conclusively that the acoustic wave interference condi¬ 
tion shown in Figure 4(c) in the main text occurs in the mechanical domain due to a cancellation between the optomechanical 
and electromechanical drive terms. 

Figure [S8] shows an optomechanical nanobeam cavity with two closely spaced mechanical modes, one the breathing mode 
which has both strong optomechanical and electromechanical coupling and a second mode which has strong electromechanical 
coupling but weak optomechanical coupling. We use the exact same procedure as in Figure 4(c) to probe the device. Starting 
with no RF power applied to the IDT, we see a single Lorentzian peak corresponding to OMIT (black curve) wherein only the 
mode with strong optomechanical coupling is visible. As we increase the RF power to the IDT keeping the phase set for acoustic 
wave interference, we see that the first mode decreases in amplitude whereas the second (higher frequency) mode rises. This 
opposite behaviour can be explained by noting that the first mechanical mode (breathing mode) is being driven by two terms, 
the optical drive term due to the beating between the carrier and the phase modulated sideband and the acoustic term due to the 
surface wave incident on the cavity. In contrast, the second mode is primarily driven by the propagating acoustic mode. Thus, 
the first mode clearly demonstrates the acoustic cancellation effect when the two driving terms are comparable in amplitude 
but opposite in phase whereas the second mode does not show any acoustic wave interference effect because there is only one 
driving term. 



FIG. S8: Normalized probe sideband transmission amplitude as a function of increasing RF power to the IDT. The lower frequency mode 
(with strong optomechanical coupling) shows a clear signature of acoustic coherent population trapping whereas the higher frequency mode 
does not. The OMIT and EMIT cases, corresponding to zero and maximum RF power applied to the IDT, respectively, are shown separately 
for clarity on the right, while the acoustic coherent population trapping case (CPT) is shown separately on the left. 
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